Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 2 April 2008 (MN MbX style file v2.2) 



IR Contamination in Galactic X-Ray Novae 



Mark T. Reynolds^*, Paul J. Callanan^, Edward L. Robinson^ and Cynthia S. Froning^ 

^Physics Department, University College Cork, Ireland 



oo 

o 

O 

<\ 



O 



> 

O 

cn 
o 

^' 

o 

oo 

o 



^Department of Astronomy, University of Texas at Austin, Austin, TX 78712 
Center for Astrophysics and Space Astronomy, University of Colorado 



2 April 2008 



ABSTRACT 

The most widely used means of measuring the mass of black holes in Galactic binaries - 
specifically the X-ray novae - involves both radial velocity measurements of the secondary 
star, and photometric measurements of its ellipsoidal variability. The latter is important in 
constraining the inclination and mass ratio, and requires as direct a measure of the flux of the 
secondary as possible. Up to now, such measurements have been preferentially carried out in 
the NIR (1 - 2.5/im), where the flux from the cooler secondary is expected to dominate over 
that from the accretion disc. However, here we present evidence of a significant non-stellar 
contribution to the NIR flux in many of those quiescent X-ray novae that are thought to contain 
a black hole primary. We discuss origins of this excess and the effect of such contamination 
on Galactic black hole mass measurements. 

Key words: stars: individual (XTE J1118H-480, GRO J0422H-32, A0620-00, GS 2000H-25, 
GS 1 124-683, GS 2023H-338, Cen X-4) - X-rays: binaries 



1 INTRODUCTION 

X-ray novae (XRNe) are a subset of the low mass X-ray binaries 
(LMXBs), containing a compact object primary (Mi ) accreting ma- 



terial via Roche lobe overflow from a low mass (M2 



1 Mq) 



, secondary star. The secondary is normally a K or M-type main se- 
quence star, whereas the primary is either a neutron star or a black 
hole with a characteristic mass of 1.4 Mq or ~ 10 Mq respec- 
tively. These systems undergo quasi-periodic outbursts (recurrence 

' timescales ~ tens of years) during which their broadband luminosi- 
ties (X-ray - radio) increase by many orders of magnitude. They 
have proven to be a fertile hunting ground for the discovery of black 
holes, with ^ 75% of the known s tellar mass black holes residing i n 
these XRNe. See the revie ws by iMcClintock & Remillard (lOOeh . 
ICharles & Coe (2006h and iFender (20061) for further information 
about the X-ray, UV/optical/IR and radio properties of these sys- 
tems respectively. 

C onclusive identification o f the primary as a black hole (Mi > 
3 MQ. lKalogera & Bavm 19961) requires further observations once 
the system has returned to quiescence. The following 3 quantities 
are required to measure the mass of the primary: (i) radial veloc- 
ity of the secondary star, K2 (ii) mass ratio, q = and (iii) 
orbital inclination, i jCasares 200^ . The first two measurements 
can be made via optical spectroscopy in quiescence. Measurement 
of the third quantity is more complicated: in most cases one must 
model the ellipsoidal modulation of the gravitationally distorted 
secondary star. Crucially, this modelling requires that one must ac- 
count for any non-stellar flux present in the lightcurve. 
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Observations of XRNe have predominately been carried out 
in the optical. In this wavelength range (3000A < A < 9OOOA), 
the accretion disc is known to contribute significantly to the ob- 
served flux. In quiescence, prominent emission lines from the disc 
are observed (i.e. IIa,II/3,HeI etc.) superposed on the spectrum 
of the secondary star. Doppler imaging of XRNe in quiescence 
shows that the a ccretion disc is still present in this state (e.g. 
iMarsh et al. 1994h . and reveals the presence of the characteristic 
hotpot where the accretion stream impacts the disc. Optical pho- 
tometry also reveals the lijrhtcurve ofjjuiescent XRNe to contain 
significant var iability e.g. iHvnes et al. (2003h : IZurita et al. (2003h : 
iNielsen (2007h . Hence, these observations show that mass trans- 
fer from the secondary star is continuously taking place. The 
contribution from the hot accretion disc is typically observed to 
decrease from the U-band (~ 3650 A) to the R-band 6500 
A). As such it has been assumed that the accretion disc should 
contribute even less at near infrared wavelengths (NIR : 1 - 
2.5^m). Initial NIR spectros copy of GS 2023+338 an d A0620- 
00 appeared to con firm this |Shahbaz et al. 199611 199 9*. although 
see Froning & Robins on 200 1|) . This appeared to justify the use of 
NIR observations, rather than optical, in determining the ellipsoidal 
variability (and constraining q & i). 

However, rece ntly it has become app arent that this might not 
always be the case. iRevnolds et al. (2007h analysed quiescent K- 
band photometry of GRO J0422+32. They found the NIR light 
curve to be dominated by a flickering component, attributed to 
emission from the accretion disc, and not an ellipsoidal modula- 
tion as one would expect if the secondary star was the dominant 
source of the observed NIR flux. Using this data, the mass of the 
black hole in this binary was estimated to be > 10.4 Mq. 
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Table 1. General properties of the various XRNe examined in this paper. See the papers indicated in column 1 for the general properties of each system. The 
veiling estimates are from the papers indicated in final column (see also Charles & Coe 2006). 



Source 


Spectral type Porb 


Distance 


Extinction 


Mx 


R-band veiling 




[hrs] 


[Kpc] 


E(B-V) 


[Mq ] 





XTEJl 118+480 " 


K5V - MOV 


4.1 


1.8 ± 0.6 


0.02 ± 0.006 


8.53 ±0.6 


0.45 


GRO J0422+32 


MOV - M2V 


5.1 


2.7 ± 0.3 


0.3 ± 0.1 


~ 10 


0.39 ' 


A0620-00 " 


K3V - K7V 


7.8 


1.2 ±0.4 


0.39 ± 0.02 


9.7 ± 0.6 


O.IO^'*^ 


GS 2000+25 


K3V - K6V 


8.3 


2.7 ± 0.7 


1.5 ± 0.1 


5.5-8.8 


0.32 


GS 1124-683 = 


K3V - K5V 


10.4 


5.5 ± 1.0 


0.30 ± 0.05 


6.95 ± 0.6 


0.15 ' 


Cen X-4 f 


K3V - K7V 


15.1 


0.9 - 1.7 


0.1 ±0.05 


1.5 ± 1.0 


0.25 ™ 


GS 2023+338 9 


KOIV 


156 


4+2.0 
^-1.2 


1.0 ±0.1 


121^ 


0.16 " 



REFERENCES: "iGehno et al. 2006t ''[Reynolds et al. 20071; '^iFroning et al. 2007t ''lloaiinou et al. 20041; "iGehno et al. 200l3; J ' iD'Avanzo et al. 2003 ; » 
IShahbaz et al. 1994 '' iTon-es et al. 2004 ' iHarlaftis et al. 19991; .' IShahbaz et al. 20041 ; '' iMarsh et al. 1994 ' ICasares et al. 1997t "" iTon-es et al. 20021; " 

ICasares et al. 19931 



Previously, iGelino et al. (20031) had measured the mass of the 
black hole in this system to be 3.97 ± 0.95 Mq; crucially they 
assumed the contribution of the accretion disc to be negligible. A 
similar result was found in an i nvestigation of the N IR flux in the 
proto-typical XRN A0620-00 bv lFroning et al. (2007h . Using mod- 
erate S/N spectra, they measure the contribution of the secondary 
star to the observed flux to be 82 ± 2% in the H-band 1.65/im). 
From this they constrain the mass of the blac k hole to be Mx = 9.7 
± 0.6 Mq. In contrast iGelino et al. (2001 lj) measured the black 
hole to have a mass of 11 ± 1.9 Mq, and in this case the secondary 
was again assumed to contribute essentially 100% of the NIR flux. 
A number of other black hole mass estimates have been made in 
the NIR (see Table [TJ and it is possible that the mass estimates for 
these systems may also be subject to additional uncertainty. 

The presence of such a contribution from a cool accre- 
tion disc component is expected t heoretically; howe ver, observa- 
tional evidence is difficult obtain jHvnes et al. 20051 ). In addition, 
at longer wavelengths (i.e. MIR: 3 - 24/xm), Spitzer observa- 
tions of a number of XRNe indicate the presence of a circumbi- 
nary disc which w ill contribute significantly at these wavelengths 
jMuno et al. 20061) . 

These results have lead us to re-examine the issue of the non- 
stellar contribution to the NIR flux of quiescent XRNe and in par- 
ticular to re-evaluate the degree to which NIR studies improve our 
chances of constraining q, i, and hence the black hole mass, in 
comparison to optical studies. These results are presented herein 
as follows: we introduce our sample systems and the relevant data 
in §2, and show the relevant spectral energy distributions (SEDs) 
in §3. It is immediately apparent that there is a significant non- 
stellar contribution to the observed flux in the NIR. Modelling the 
multi-wavelength SEDs shows that this excess is consistent with 
blackbody emission from a combination of the accretion disc and a 
circumbinary disc. A discussion of our results follows in §4. 



2 DATA 

Our sample consists of XRNe with quiescent data spanning the op- 
tical to at least the NIR wavelength range. In addition, we choose 
only those systems that contain main sequence companions (class 
V), with the exception of GS 2023+338. This XRN has a much 
longer orbital period and a sub-giant companion star (class IV). 
However, as this is the only other quiescent XRN with cover- 
age at Spitzer wavelengths, it has been included for completeness. 



The XRNe that make up our sample include 6 suspected to con- 
tain a black hole primary and a single neutron star system. In Ta- 
ble [T] we list the primary system parameters for each binary. Our 
data consists of archival published optical and NIR photometry as 
well as a number of mid- IR (MIR: 3 - 24 pm) data points from 
the Spitzer observatory ( Muno et al. 20061) . In particular Spitzer 
photometry exists for 4 of the systems, namely: A0620-00, XTE 
Jl 118+480, GS 2023+338 & Cen X-4. We have also recently ob- 
tained a medium resolution K-band spectrum of the black hole sys- 
tem GS 2000+25, which we discuss below. In Table[2l we list the 
dereddened flux from each system in the wavelength regions of in- 
terest. 

2.1 GS 2000+25 Spectroscopy 

K-band spectra of GS 2000+25 we re obtained with the Near- 
IR imager (NIRI, Hodapp et al. 2003) using the f/6 camera at 
the Gemini-North telescope. Medium-resolution spectra of GS 
2000+25 were obtained on the nights of 2006 August 23, 25 and 
September 1 UT as part of the Gemini service program (Prog, id: 
GN-2006B-Q-85, PI: Robinson). The K grism (1.9 - 2.49/im) was 
used in conjunction with a 0.23 " slit, giving a resolution R ~ 1 100. 
Individual exposure times were 300s, dithered to 5 positions along 
the slit. In total 37 spectra of the system were acquired yielding a 
total exposure time of ~ 3hrs. At the time of our observations pho- 
tometry revealed the K-band magnitude of GS 2000+25 to be 16.6 
± 0.15, where the error is dominated by the error in the 2MAS^| 
stars relative to which our image is calibrated. This is consistent 
with previous measurements of the system, showing it to be in qui- 
escence at the time of our observations. 

The spectra were observed in 3 separate groups of 23, 8 
and 6 science frames taken on each of the three nights outlined 
above. The airmass was typically in the range 0.95 - 1.15 dur- 
ing the observations as such slitlosses due to atmospheric differen- 
tial refraction were negligible. The one-dimensional spectra were 
extracted within IRAI0 using both the NOAO.ONEDSPEC and the 

^ The Two Micron All Sky Survey is a joint project of the University of 
Massachusetts and the Infrared Processing and Analysis Center, California 
Institute of Technology, funded by NASA and the National Science Foun- 
dation. 

^ IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
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2x10^ 2.1x10* 2.2x10* 2.3x10-' 2.4x10-' 

Wavelength [ A ] 

Figure 1. Gemini NIRI K-band spectmrn of GS 2000+25 (top), the Hel and 
Br7 emission lines are indicated, also indicated are the positions of promi- 
nent atornic absorption lines and the CO-bandheads. The crosses mark the 
position of telluric residuals, the spectrum of a K5V star is plotted under- 
neath for comparison. 



any CO bandheads in the GS 2000-1-25 spectrum. As the S/N de- 
grades further at longer wavelengths, our efforts were confined to 
the CO (2-0) band head at 2.294/im. We used a model of this band- 
head in the K5V star spectrum to place an upper limit on its pres- 
ence in the observed GS 2000-1-25 spectrum. We can rule out the 
presence of this bandhead with an equivalent width of > 20% of 
that expected from a K5V star (solar abundance) at the 99% confi- 
dence level. 

The spectral slope is observed to be consistent with that ex- 
pected from the K5V secondary star in this system for A ^ 2.3/im: 
for longer wavelengths, the slope appears to flatten. In an effort to 
constrain the amount of non-stellar flux present a number of spec- 
tra were simulated. Noise was added to a K5V spectrum consistent 
with the secondary in this system until a S/N ratio similar to the 
science spectrum was achieved. To this a constant contribution rep- 
resenting an accretion disc with a temperature profile of the form T 
oc i"^''^ was added (see We find for a fractional disc contri- 
bution of 25%, the simulated spectral slope is no longer consistent 
with the observed spectrum, indicating a likely disc contribution of 
< 25% in the 2.0 - 2.3/im wavelength region and more at wave- 
lengths greater than this. 



GEMINI. NIRI packages. Both sets of extracted spectra were found 
to be consistent. Unfortunately, due to the nature of the service ob- 
servations a significant fraction of the spectra were taken in poor 
to deteriorating conditions resulting in a large number of unusable 
exposures. Hence the final signal to noise ratio (S/N) of our sci- 
ence spectrum was lower than anticipated. Spectra of the AOV stars 
HD182761 and HD197291 were also taken before and after the 
science frames respectively to aid with the removal of the telluric 
features and flux calibration of our science spectrum. 

Correction of the science frames fo r telluric absorptio n was 
carried out following the prescription of IVacca et al. (2003h . The 
low S/N of the individual spectra complicated the removal of the 
telluric features, with a number of residuals remaining. The re- 
sulting spectra were corrected for the Doppler shift of the spec- 
tral lines induce d by the orbital motio n of the secondary star using 
the ephemeris of lloannou et al. (20041) and median combined so as 
to negate the effect of cosmic-ray hits and other spurious effects, 
while at the same time maximising the available S/N. This spec- 
trum was then dereddened using the extinction listed in Table[T] As 
the slit losses are large, only a relative flux calibration was possi- 
ble. The resulting spectrum (~ 2hrs exposure) is displayed in Fig. 
[TJ emission lines of Hel and Br7 are observed. 



3 ANALYSIS 

3.1 GS 2000-K25 K-band Spectrum 

The K-band spectrum in Fig. [T] is displayed at a resolution of ~ 
1100 (270 km s"^) and has a S/N ratio of ~ 10 at 2.3/im. We 
observe emission features consistent with both Hel and Br7 with 
equivalent widths of ~ 6 A & 12A respectively. Also indicated are 
the positions of a number of the atomic absorption features ex- 
pected from the secondary star. 

We also attempted to place on upper limit on the presence of 



tronomy Inc., under cooperative agreement with the National Science Foun- 
dation. 



3.2 Multiwavelength Spectral Energy Distributions 

In Figures |2] & [S] we plot the extinction corrected optical/IR 
SEDs of the XRNe as listed in Table |2l where the error bars 
account for the uncertainty in the photometry and the redden- 
ing. Here, we fit model atmospheres of the appropriate spec- 
tral type to the SEDs of the quiescent XRNe: these data in- 
clude optical, NIR and Spitzer measurements. In contrast to previ- 
ous authors who normalised to the H-band (Gelino et al. 2006) or 
K-band flux dMuno et al. 20061 : iGelino et al. 2003ll 2001^ , where 
the IR contamination from the accretion disc was assumed to 
be minimal, we choose instead to normalise these models rel- 
ative to the measured secondary contribution in the R-band, 
which has generally been more accurately determined using op- 
tical spectr oscopy. The models we fit are NextGen mode l at- 
mospheres taauschildt et al. 1999^ : iHauschildt et al. 1999bh cor- 
responding to the spectral type of the secondary star in each system. 
We restrict ourselves to those veiling measurements for which the 
simultaneously measured R-band flux is consistent with the value 
used in the SED. The veiling in all cases is defined as a percent- 
age of the total system flux emitted by the accretion disc, typically 
measured in the region of the Ha line (6562A): 

-r-\ Fdisc 

star ^ disc 

and is indicated in the top left hand corner of each SED in Fig. |2] 
&[3]and in the final column of Table[T] As the various photometric 
data points used were non-simultaneous, we do not attempt to rig- 
orously quantify the excess in the NIR. The quiescent magnitudes 
of XRNe are typi cally observed to vary by ~ 0. 1 magnitudes (e.g. 
IZurita et al. 20031) , in addition to the ^ 0.2 magnitude periodic vari- 
ability due to the intrinsic ellipsoidal modulation of the secondary 
star. However, when normalised this way, the SEDs of these XRNe 
are all seen to exhibit an excess of flux at NIR wavelengths compa- 
rable to that measured in the optical (also indicated in Figures |2]& 
[3]l. In this sense the NIR appears to be just as contaminated as the 
optical, in the context of measuring the flux from the secondary. 
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Table 2. The spectral energy distribution data. All fluxes are units of erg s~^ cm~^. 



Source 


B 


V 


R 


I 


J 


H 


K 


4.5^m 


8/im 


24/^r?i 


Radio 




4400A 


5500A 


6400A 


7900A 


1.25/im 


1.65/im 


2.2fim 








8.5GHz 


XTEJ1118+480 


2.54e-13 


2.89e-13 


2.94e-13 




3.50e-13 


2.78e-13 


2.1e-13 


3.07e-14 


1.69e-14 


> 2e-15 




GRO J0422+32 ^ 




7.35e-14 


1.18e-13 


1.39e-13 


2.32e-13 


2.07e-13 


1.07e-13 










A0620-00 4>5.6.7.2 


1.23e-12 


2.09e-12 


2.30e-12 


2.38e-12 


3.25e-12 


2.88e-12 


1.58e-12 


2.99e-13 


9.34e-14 


6.75e-15 


4.33e-18 


GS 2000+25 






8.34e-13 


6.34e-13 


7.70e-13 


7.52e-13 


3.53e-13 










GS 1124-683 


8.98e-14 


2.20e-13 


2.22e-13 


2.49e-13 


2.90e-13 


2.98e-13 


1.43e-13 










Cen X-4 "'^ 


6.15e-13 


9.09e-13 


1.39e-12 


1.59e-12 


2.08e-12 


1.84e-12 


l.lle-12 


1.33e-13 


3.56e-14 


> 3.8e-15 




GS 2023+338 12,13,2 


6.88e-12 


1.54e-ll 


2.34e-ll 




3.11e-ll 


2.57e-ll 


1.29e-ll 


2.01e-12 


5.44e-13 


1.91e-14 


~4e-17 



REF ERENCES: i(Gel ino et al. 2006; 2 Muno et al. 2006: '^iGelino et al. 20031 : "^iGelino et al. 2001 bl: %ronins & Robinson 2001 ; ^iHasweU et al. 199i 
lOallo et al. 2006 : ^ .Callanan & Charles 1991, ; '' .Callanan et al 19961 : l°lGehno et al. 2001al : ^i lstiahbaz et al. 199^; ^2 ,Casares et al. 19931 ; " 
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3.3 A Simple Model 

The current paradigm for understanding the emission from qui- 
escent XRNe involves a thin disc which transitions to a quasi 
spherical inner flow at a distance of ~ 10^ - 10'* Schwarzschild 
radii from the compact object. The inner flow is t hought to con- 
sist of an advection dominated accretion flow (see iNaravan 20081 
and references therein), although it has also b een argued that the 
inner region in stead consists of a jet/outflow jpender et al. 20031 : 
lYuan & Cui 200 5). In an effort to place a constraint on the source of 
the observed excess NIR flux, we have attempted to model the ob- 
served SED. The spectral energy distributions are fit with a multi- 
component model consisting of: 

(i) A spherical blackbody representing the emission from the 
secondary star. 

(ii) An accretion disc. 

The flux from the accretion disc may be modelled in the form of 
a multi-colour blackbody where the temperature profile across the 
surface of the optica lly thick geometrically thin accretion disc is 
jMitsuda et al. 1984) : 



Teff(r) 



/ SGMirii 



1/4 



(1) 



This results in Teff (r) oc r~''^*for r ^ Rin in the case of a vis- 
cously heated steady state disc. 

Alternatively, we might expect a flatter temperature profile: 
for example, Tefr(r) oc r~^''2 has been observed in the 
quiescent accretion di s c in some c ataclysmic variables (CVs, 
IWood et al. 198911 19921: [iVIenou 2002h . Such a temperature profile 
is intermediate between the classic steady state viscously heated 
accretion disc (r^"'^^) and the irradiated disc case (r"'^ *''). 

(iii) A circumbinary disc. 
If we assume the circumbinary disc to be opaque, flat and p assively 
illuminated by the central star JChiang & Goldreich 1997h . we ex- 
pect the following temperature profile: 



T.«(r)=(3^ 



1/4 



3/4 



(2) 



where R* & T* are the radius and temperature of the irradiating 
star respectively. 



3.4 Results 

Given the non-simultaneous nature of the observations we do not 
try to rigorously fit the data (e.g. via minimisation). None the 
less, in each case it is clear that a model consisting of an accre- 
tion disc extending from 10"^ - 10^ Schwarzschild radii (Rg) to an 
outer radius of ~ 0.5 the Roche lobe radius of the compact object 
(Rli), in addition to the secondary star, is required to provide the 
observed NIR flux. This is consist ent with previous observations 
of the accretio n disc in quiescenc e jMcClintock et al. 2003111993: 
iNaravan 20021 : iMarsh et al. 19941) . At Spitzer wavelengths, the 
flux appears to be a combination of emission from the accre- 
tion disc/circumbinary disc. The SEDs of all the XRNe are seen 
to follow this general pattern. We emphasise that although the 
following analysis is qualitatively similar to that carried out by 
iMuno et al. (2006b , it differs in the important aspect that we nor- 
malise to the R-band flux, which results in a significant non-stellar 
component to the observed NIR flux (see We discuss the fits 
to each system in more detail below. 



3.4.1 A0620-00 

In Figure |4] we display the resulting fit in the case of A0620-00. 
A steady state disc fit is displayed on the left: a disc with a flatter 
temperature profile (T oc r""'^) is displayed on the right. We see 
that the disc with the flatter temperature profile provides a better 
fit to the observed SED. The basic model consists of a blackbody 
at the temperature of the model atmosphere displayed in Figure 
|2l a 'flat' accretion disc extending from ~ lO** Rs to ~ 0.5 Rli 
and a circumbinar y disc consistent with the previous estimates of 
iMuno et al. (20061) . Such an accretion disc will contribute ~ 40% 
of the K-band flux. In comparison, iFroning et al. (20071) recently 
determined the accretion disc to be contributing ~ 23% of the K- 
band flux through moderate S/N spectroscopy. Given the uncertain- 
ties discussed elsewhere (§3.2, §4), we regard these values to be 
consistent with each other. 



3.4.2 XTEJl 118+480 

Here we find a similar model to that described above, in particular 
a flat temperature profile accretion disc is required. The slope of 
the SED is consistent with emission from the secondary plus the 
accretion disc out to 4.5/xm, however the 4.5 - Sfim slope clearly 
indicates the presence of an additional source of flux. 
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R-band veiling ~ 10% 
K-band veiling - 40% 
K5V secondaiy 



45|lm 
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10'= So" 



logju Frequency [ Hz ] 
(a) A0620-00 




logjp Frequency [ Hz ] 
(b) XTE J1118+480 



24^m 



K H J 



1 R V B 



R-band veiling ~ 32% 
K-band veiling ~ 25% 
K5V secondaiy 




10'* 

logjg Frequency [ Hz 
(c) GS 2000+25 



logiQ Frequency [ Hz ] 
(d) GRO J0422+32 



Figure 2. Opdcal/IR SEDs for the XRNe outlined in the text. The grey curve is the absorption connected SED, where the enor bars account for the uncertainty 
in the photometry and reddening. The black curve represents the blackbody corresponding to the temperature of the model atmosphere (Ught-grey). The R-band 
accretion disc contamination is listed (see Table [T) along with the resulting K-band contamination and the adopted spectral type of the secondary star in the 
top left. 



In contrast to lMuno et al. (200^ . we find a standard circumbi- 
nary disc is unable to account for the observed excess. This is 
because we account for the accretion disc veiling, which results 
in a lower contribution from the secondary star at Spitzer wave- 
lengths in this work, requiring a greater flux contribution from the 
circumbinary disc. Instead, we find a circumbinary disc with a flat- 
ter temperature profile similar to the accretion disc is required. In 
this case the accretion disc will contribute ~ 60% of the flux at NIR 
wavelengths. 



3.4.3 GS 2023+338 

As is immediately obvious from the SED in Fig. [3] the flux from 
this system is dominated by the emission from the type IV sec- 
ondary star. We find the accretion disc to contribute in the region of 
25% of the flux in the JHK region of the spectrum. The excess at 24 
microns is also consistent with an origin in the accretion disc. How- 
ever, we are unable to rule out a contri bution from a circumbinary 
disc as suggested bv lMuno et al. (200d) . 



3.4.4 Cen X-4 

Cen X-4 is the only system in our sample where the compact object 
is a neutron star. However, as in the other systems the secondary 
star is unable to account for the NIR/optical flux, with there be- 
ing a large excess over that from the secondary at NIR and Spitzer 
wavelengths. The magnitude of the contribution from the accretion 
disc is consistent with that observed in the black hole systems men- 
tioned above. The disc is required to extend from an inner radius of 
0.3 Rq to an outer radius of 0.5 Rli- We also note that the Spm 
flux in conjunction with the 24pm upper limit allows for the pos- 
sible presence of a circumbinary disc, which c ontributes up to 20% 
of the Sfim flux, in contrast to the analysis of iMuno et al. (200^ . 
Such a circumbinary disc will not contribute to the observed K- 
band emission. 



3.4.5 Systems Without Spitzer Measurements 

GRO J0422+32: Previous observations of GRO J0422-I-32 high- 
lighted the likely presence of a large non-stellar NIR component in 
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Figure 3. As in Figure|2] 



this system jRevnolds et al. 2007h . Modelling of the SED supports 
this claim, with the accretion disc component providing ~ 40% of 
the observed NIR flux. 

GS 2000+25: Again, here we find the accretion disc contribution 
to be large, ~ 30% in the K-band. 

GS 1124-683: The flux from the accretion disc is seen to be ~ 40- 
50% of the total flux at NIR wavelengths. 



4 DISCUSSION 

As has been demonstrated in the previous section, it is clear that the 
NIR flux observed from quiescent X-ray novae contains a signifi- 
cant component that does not originate in the secondary star. From 
an analysis of the spectral energy distributions of these XRNe, it 
now appears that this NIR excess is dominated by thermal emis- 
sion from the cooler outer regions of the accretion disc. This excess 
requires the presence of a cool accretion disc with a temperature 
profile T (X r^"'^; such a cool disc comp onent (T ~ 3000 - 4000 
K) is expected on theoretical grounds (e.g. lHvnes et al. 2005h . This 
temperature profile is intermediate between the classic steady state 
viscously heated accretion disc (r^"'^^) and the irradiated disc case 
(r-o «). The use of a temperature profile of this form is supported 



by the observation of similarly shallow temperature profiles from 
ecli pse mapping of the accretion discs in quiescent dwarf novae 
(see lMenou 20021 and references therein). 

Given the difficulty in accurately determining the spectral type 
of the secondary star in XRNe, we also investigated the effect of 
varying the spectral type (within the range listed in TablefTJ on the 
magnitude of the NIR excess. In particular, a later spectral type 
could account for a significant portion of the observed NIR ex- 
cess (Fig.[2]&[3](. Using the latest secondary star spectral types al- 
lowed, as previously determined (see Table[TJ, the analysis in H3.4\ 
was repeated. The observed NIR flux is found to remain system- 
atically greater than the relevant model atmosphere for all the sys- 
tems in our sample, although the magnitude of the excess is found 
to decrease as expected. In the 7 systems in our sample, the ob- 
served K-band excess will remain high ranging from ^ 10% in 
A0620-00 (K7V secondary) to 40% in XTE Jl 1 18-1-480 (MOV sec- 
ondary). The veiling in each system is typically known to within ± 
5%, which can lead to an error in the disc flux ranging from 50% 
(A0620-00: 10 ± 5) to as little as 12% (Cen X-4: 25 ± 3). We adopt 
this as an appropriate estimate for the accuracy of our K-band veil- 
ing measurements. Regardless of the magnitude of the excess, it is 
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the systematic presence of this excess in all of the surveyed systems 
that is the most convincing evidence for its existence. 

In the following sections, we discuss additional spectroscopic 
evidence for accretion disc contamination of the NIR flux in XRNe, 
which support the SED analysis given in i]3.4l and above. In addi- 
tion, we consider possible jet emission in the NIR, briefly discuss 
the MIR flux and compare the systems containing black hole pri- 
maries to the neutron star system. We end by considering the im- 
plications of neglecting to properly account for this additional flux 
when determining black hole masses. 



4.1 Previous K-band Spectroscopy 

NIR spectroscopy of quiescent black hole XRNe has also pro- 
vided evidence for a non-negligible contribution from the accre- 
tion disc. Recent mo derate resolution observations of A0620-00 by 
iFroning et al. (2007h . support the presence of a near-infrared excess 
in this system. Here the spectral slope was observed to deviate sig- 
nificantly from that expected of the secondary star alone. Detailed 
analysis showed the secondary star to cont ribute 82 ± 2% of the 
H-band flux and ~ 77% of the K-band flux. [Reynolds et al. (200"7h 
analysed low resolution spectra of the XRN GRO J0422-I-32 that 
displayed significant Br7 emission; in combination with photomet- 
ric observations they showed the accretion disc to be contributing 
up to ~ 30% or more of the observed K-band flux. 

4.2 The GS 2000+25 Spectrum 

In this work, we have obtained a medium resolution K-band spec- 
trum of GS 2000+25. Emission lines consistent with the accretion 
disc are detected (Hel, Br7), whereas the intrinsic absorption lines 
of the secondary star appear to be absent (although the S/N is low, 
see Fig.[T]l. Modelling of the CO (2-0) bandhead at 2.29Apm al- 
lows us to place a limit on the equivalent width of this feature in 
our spectrum. We can exclude the presence of this bandhead with 
an equivalent width of ^ 20% of that expected from a K5V star, at 
the 99% confidence level. On the other hand, analysis of the spec- 
tral slope of our spectrum allows us to conservatively constrain the 
accretion disc to be contributing < 25% of the flux for A ^ 2.3/im, 
and more at wavelengths greater than this, although a more accu- 
rate determination is limited by our low S/N ratio. We regard this 
as consistent with the value estimated in i]3.2l & Fig.|2] 

These results indicate that the CO abundance in the secondary 
star in GS 2000+25 is anomalous. If the CO abundance were sim- 
ilar to solar, we would expect to detect the CO (2-0) bandhead at 
approximately 75% of the equivalent width found in the K5V spec- 
trum, with the remainder of the line being filled in by the accretion 
disc continuum, which we have shown contributes < 25% of the 
observed flux. Instead, we limit the maximum CO (2-0) bandhead 
strength to be 20% of that expected from a K5V star at solar abun- 
dance. 

In recent observations of A 0620-00 at similar res- 
olution (albeit hig her S/N) both iFroning et al. (2007h and 
iHarrison et al. (20"07b detect the CO bandheads from the secondary 
star in absorption. The CO lines detected in A0620-00 are observed 
to be anomalous ly weak, hinting at possible CNO processing in the 
secondary star jFroning et al. 2 007). Observations of the XRNe 
XTE Jl 1 18+480 in outburst reveal a similar situation, with UV 
spectra displaying evidence for anomalous abundances of C, N 
& O. This p oints to CNO proce ssing having also taken place in 
this system jHaswell et al. 2002 b . As explained above, if the C 



abundance in the secondary star in GS 2000+25 was anomalous 
as is the case in A0620-00 ([C/H] = -1.5), this would provide an 
explanation for their non-detection. 

We also not e the K-band spectroscopy of WZ Sge 
jHoweU et al. 2004h in which H2 and CO-bandhead emission was 
detected. These lines are thought to be produced in dense cool re- 
gions as o ne might expect in the outer regions of the accretion disc 
in XRNe ( IHvnes et al . 2005). Hence, it is possible that our non- 
detection of the secondary star CO-bandheads in absorption is due 
to a combination of an anomalous carbon abundance and infilling 
of the CO-bandheads by similar emission lines emanating from the 
outer regions of the accretion disc. Unfortunately definitive con- 
clusions require high er S/N observa tions as demonstrated in case 
of A0620-00 bv .Froning et al. (2007h . 

4.3 Non-Thermal Jet Emission 

In an effort to constrain any possible jet contribution to the ob- 
served SEDs, a jet component of the form oc was added to 
the model in |3.4| The jet is assumed to be flat and optically thick 
(a = 0), extending from radio wavelengths and breaking to the op- 
tically thin re gime (a < 0) in the optical/IR wavelength range (i.e. 
iFender 200^ . 

There does not appear to be evidence for a significant non- 
thermal contribution at NIR wavelengths in any of the systems in 
our sample. In A0620-00 and GS2023+338, the jet is seen to con- 
tribute at most a few percent of the observed fl ux. Even given the 
uncertainty in th e Spitzer detections at 2Ap,m (i Muno et al. 20061 : 
iGallo et al. 20071) . a jet will contribute ^ 5% and < 10% of the K- 
band flux in A0620-00 and GS 2023+338 respectively. It appears 
that jet emission in quiescence only becomes an appreciable per- 
centage of the emitted fl ux as we proceed t o the mid-IR region of 
the spectrum as noted bv lMuno et al. (200^ . 

4.4 MIR emission 

As we proceed to longer wavelengths, the observed flux will be- 
come dominated by emission from the accretion disc/circumbinary 
disc with a possible contribution also from a jet, in agreement 
with the analysis of Muno et al. (2006). In the case of the sources 
with mid-IR spectral coverage (A0620-00, XTE Jl 118+480, GS 
2023+338, Cen X-4), the Spitzer data favours the presence of a 
circumbinary disc contribution to the observed flux. Specifically, 
whereas a circumbinary disc is capable of accounting for all the 
mid-IR flux, a flat jet alone cannot. 

4.5 Black Hole vs Neutron Star XRNe 

Observations of quiescent XRNe with the Chandra X-ray observa- 
tory have discovered that the X-ray luminosity of systems contain- 
ing black hole primaries appear to be syst ematically fainter than 
those containing neutron stars (see I Kong et al. 200Z : iLasota 20071 
and references therein). This is thought to be due to the fact that 
the accreting matter strikes the hard surface in the case of the neu- 
tron star whereas if the primary is a black hole, the matter is sim- 
ply advected across the event horizon (although an alternate in- 
terpretation envisions the matter being expelled in a jet/outflow, 

t where the flux is measured in units of milli-Janskys: Imjy = 10^'^'^ erg 
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Figure 4. Blackbody fit to the A0620-00 SED consisting of (i) star - solid line (ii) multi-colour blackbody disc - dashed line (iii) circumbinary-disc - dotted 
line. The observed SED is in grey with the model SED indicated by the thick line. 



iFender et al. 2003h . Hence, it is worth asking if there is any ap- 
preciable difference in the SEDs of the sources in our sample. 

Even though we only model a single system containing a neu- 
tron star primary, there does not appear to be any discernable dif- 
ference between the spectral energy distribution of Cen X-4 and the 
6 systems containing a black hole primary. This is what one would 
expect if the emission in the optical/NIR region is dominated by the 
secondary star and the cooler outer regions of the accretion disc in- 
stead of radiation from processes taking place closer to the compact 
object. 

4.6 Black Hole Mass Estimates 

In light of the fact that the NIR flux from XRNe appears to contain 
a significant non-stellar contribution, it is worth reconsidering pre- 
vious black hole mass estimates made via IR observations, which 
neglected to include this non-stellar flux in the analysis. We can es- 
timate the effect of this additional co mponent using the Eclipsing 
Light Curve code (ELc lOrosz et al.^OOOi) in combination with the 
previously determined system parameters. We then add a contribu- 
tion from an accretion disc with a temperature profile of the form 
T(r) oc r^"'^ and the inner and outer disc radii set so as to agree 
with the values determined in i]3.4| We estimate that for an accre- 
tion disc contribution of only 20%, the black hole masses measured 
via NIR photometry have generally been overestimated by 1 - 2 
Mq in each case; this increases to as much as 4 Mq in the case of 
a 50% disc contribution. 



5 CONCLUSIONS 

We have shown that there is a significant amount of contamina- 
tion present from the accretion disc/circumbinary disc in the IR 
portion of the SED of quiescent XRNe. Fits to the SEDs reveal 
NIR excesses in each of the 7 systems studied. We have also pre- 
sented new K-band spectroscopy of GS 2000-1-25, which also shows 
some evidence for contamination: this joins two other black hole 
XRNe with confirm ed IR excesses from s pectroscopy (A0620-0 - 
iFroning et al. 20071 and ORG J0422-I-32 - iRevnolds et al. 2007h . 



Based on these results, we believe that the currently accepted 
paradigm, in which the ellipsoidal variations at NIR wavelengths 
are assumed to be undiluted by other sources of flux in the binary, 
is not valid. We conclude that the NIR offers no significant advan- 
tage over optical observations in the measurement of ellipsoidal 
variability (and the determination of mass ratio and orbital incli- 
nation). Indeed, assuming the contrary introduces large, systematic 
errors in the mass estimates for the compact objects in these binary 
systems. 
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